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ABSTRACT 

Several gamma-ray burst (GRB) luminosity indicators have been proposed, which can 
be generally written in the form of L = c where c is the coefficient, Xi is the j-th 

observable, and a* is its corresponding power-law index. Unlike Type-la supernovae, 
calibration of GRB luminosity indicators using a low-redshift sample is difficult. This 
is because the GRB rate drops rapidly at low redshifts, and some nearby GRBs may 
be different from their cosmological brethren. Calibrating the standard candles using 
GRBs in a narrow redshift range (Az) near a fiducial redshift has been proposed re¬ 
cently. Here we elaborate such a possibility and propose to calibrate {cii} based on the 
Bayesian theory and to marginalize the c value over a reasonable range of cosmological 
parameters. We take our newly discovered multi-variable GRB luminosity indicator, 
E- lso = cEp 1 t , £ 2 , as an example and test the validity of this approach through simula¬ 
tions. We show that while c strongly depends on the cosmological parameters, neither 
ai nor 02 does as long as A z is small enough. The selection of A z for a particular GRB 
sample could be judged according to the size and the observational uncertainty of the 
sample. There is no preferable redshift to perform the calibration of the indices {ai}, 
while a lower redshift is preferable for c-marginalization. The best strategy would be 
to collect GRBs within a narrow redshift bin around a fiducial intermediate redshift 
(e.g. 2 C ~ 1 or z c ~ 2), since the observed GRB redshift distribution is found to peak 
around this range. Our simulation suggests that with the current observational preci¬ 
sion of measuring GRB isotropic energy (U; so ), spectral break energy ( E p ), and the 
optical temporal break time (£&), 25 GRBs within a redshift bin of A z ~ 0.30 would 
give fine calibration to the Liang-Zhang luminosity indicator. 
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1 INTRODUCTION 

The cosmological nature (Metzger et al. 1997) of long 
gamma-ray bursts (GRBs) and their association with star 
formation (e.g., Totani 1997; Paczynski 1998; Bromm & 
Loeb 2002) make GRBs a new probe of cosmology and 
galaxy evolution (e.g. Djorgovski et al. 2003). Gamma-ray 
photons with energy from tens of keV to MeV from GRBs 
are almost immune to dust extinction, and should be de¬ 
tectable out to a very high redshift (Lamb & Reichart 2000; 
Ciardi & Loeb 2000; Gou et al. 2004; Lin et al. 2004). Sev¬ 
eral plausible GRB luminosity indicators have been pro¬ 
posed, including luminosity-variability relation (Fenimore 
& Ramirez-Ruiz 2000; Reichart et al. 2001), luminosity- 
spectral lag relation (Norris, Marani, & Bonnell 2000), stan¬ 
dard gamma-ray jet energy {E 1 j e t ){Frail et al. 2001; Bloom 
et al. 2003), isotropic gamma-ray energy ( Ei so ) - peak spec¬ 
tral energy ( E p ) relation (Amati et al. 2002), E lt j e t — E v 


relation (Ghirlanda et al. 2004a), and a multi-variable re¬ 
lation among E lso , E p and the break time of the optical 
afterglow light curves ( tb ) (Liang & Zhang 2005). Attempts 
to use these luminosity indicators to constrain cosmologi¬ 
cal parameters have been made (e.g. Schaefer 2003; Bloom 
et al. 2003; Dai, Liang, & Xu et al. 2004; Ghirlanda et al. 
2004; Friedmann & Bloom 2005; Firmani et al. 2005; Liang 
& Zhang 2005; Xu, Dai, & Liang 2005; Xu 2005; Mortsell 
& Sollcrman 2005; Wang & Dai 2006). With the discovery 
of the tight Ghirlanda-relation (Ghirlanda et al. 2004a) and 
the more empirical LZ-relation (Liang & Zhang 2005), it is 
now highly expected that GRBs may become a promising 
standard candle to extend the traditional Type-la SN stan¬ 
dard candle to higher redshifts (e.g. Lamb et al. 2005). 

In order to achieve a cosmology-independent standard 
candle, one needs to calibrate any luminosity indicator. Oth¬ 
erwise, one inevitably encounters the so-called “circularity 
problem” (e.g. Firmani et al. 2004; Xu et al. 2005 for dis- 
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cussion). In the case of Supernova cosmology, calibration is 
carried out with a sample of Type-la SNe at very low red- 
shift so that the brightnesses of the SNa are essentially in¬ 
dependent on the cosmology parameters (e.g., Phillips 1993; 
Riess et al. 1995). In the case of GRBs, however, this is very 
difficult. The observed long-GRB rate falls off rapidly at 
low redshifts, as is expected if long-GRBs follow global star 
formation. Furthermore, some nearby GRBs may be intrin¬ 
sically different. Observations of GRB 980425, GRB 031203 
and some other nearby GRBs indicate that they differ from 
typical GRBs by showing low isotropic energy, simple light 
curve, large spectral lag and dimmer afterglow flux (e.g., 
Norris 2002; Soderberg et al. 2004; Guetta et al. 2004; Liang 
& Zhang 2006). Although some GRBs in the optically-dim 
sample of Liang & Zhang (2006) still follow the Ghirlanda- 
relation and LZ- relation, some very nearby events (e.g. GRB 
980425 and GRB 031203) are clearly outliers of the Amati- 
relation and also likely the outliers of the Ghirlanda-relation 
and the LZ-relation. If at least some low-redshift GRBs are 
different from their cosmological brethren, it is very difficult 
to calibrate the GRB standard candle using a low-redshift 
sample. 

Recently the possibility of calibrating the standard can¬ 
dles using GRBs in a narrow redshift range (Az) near a fidu¬ 
cial redshift has been proposed (Lamb et al. 2005; Ghirlanda 
et al. 2005) 1 . In this paper we elaborate this method (Lamb 
et al. 2005; Ghirlanda et al. 2005) further based on the 
Bayesian theory. We propose a detailed procedure to cal¬ 
ibrate {di} with a sample of GRBs in a narrow redshift 
range (Az) without introducing a low-redshift GRB sample 
and marginalize the c value over a reasonable range of cos¬ 
mological parameters. The method is described in §2. We 
take our newly discovered GRB luminosity indicator as an 
example to test the approach through simulations (§3). The 
results are summarized in §4 with some discussion. 


2 CALIBRATION METHOD 

A GRB luminosity indicator can be generally written in the 
form of 

L{Q)=c(Cl)Q(Cl\X;A), (1) 

where c(fi) is the coefficient, 12 is a set of cosmological 
parameters, and Q(Cl\X-,A) is a model of the observables 
X = {xi} (measured in the cosmological proper rest frame) 
with the parameter set A = {a;}, which is generally written 
in the form of Q(Cl\X-, A) = Y[ x i'\ n - Since in GRB lumi¬ 
nosity indicators the parameters { Xi } are usually direct ob¬ 
servables (e.g. E p , tb, etc) that only depend on z but not on 
the cosmological parameters, the above expression naturally 
separates the Sl-dependent part, c(fi), from the fi-insensitive 
part, Q(Cl\X- A). This allows us to develop an approach to 
partially calibrate the luminosity indicators without requir¬ 
ing a low-redshift GRB sample. Our approach is based on 
the Bayesian theory, which is a method of predicting the 
future based on what one knows about the past. Our cali¬ 
bration process can be described as follows. 


1 The similar idea was also discussed in an earlier version of Liang 
&; Zhang (2005). 


(1) Calibrate A using a sample of GRBs that satisfy 
a luminosity indicator and are distributed in a narrow red¬ 
shift range zo € z c ± Az. Luminosity distance as a function 
of redshift is non-linear, and the dependence of the luminos¬ 
ity distance on the cosmology model at different redshift is 
different. Such a sample reduces this non-linear effect. The 
parameter set A can be then derived by using a multiple 
regression method in a given cosmology Cl, A(Cl,zo). The 
goodness of the regression is measured by Xmm(^> z o ), 

2 ^ [logL ! (»,zo)-logL ! (n,z 0 )] 2 

Xminw", Zo) 2 ' (2) 

i log L l (Cl,z 0 ) 

where N is the size of the sample, <x log ii^ ZQ s ) is the error 
of the empirical luminosity from the observational errors of 
observables, and logL l (f2,zo) is the theoretical luminosity. 
The smaller the reduced Xmin(^>zo)i the better the regres¬ 
sion, and hence, the higher the probability that A(f2, zo) is 
intrinsic. Assuming that the Xmm(^X-Zo) follows a normal 
distribution, the probability can be calculated by 

P(Cl,zo) oc e _Xmin ^’ z °^ 2 . (3) 


The calibrated A with a sample distributed around zo is 
thus given by 

In z o)P(&, zo)dO. 

/ n P(n,z 0 )dn ’ () 

and its root mean square ( rms ) could be estimated by 


^ / n [A(H, zo) - A(z 0 )] 2 P(Cl, z 0 )dCl 

f n P(Q,zo)dn 


( 5 ) 


where Ao(zo) is the unweighted mean of A(fi, zo) in different 

Cl. 

(2) Marginalize the c value over a reasonable range for 
a given GRB sample. The c value depends strongly on the 
cosmological parameters, so it can only be calibrated with a 
low redshift sample. Because of the reasons discussed above, 
such a low-z sample is hard to collect. We therefore do not 
calibrate the c value but rather marginalize it over a rea¬ 
sonable range of cosmological parameters for a given GRB 
sample. For a given c one can derive an empirical lumi¬ 
nosity L{Cl, c, Ao, zo) from the luminosity indicator and its 
error. The x 2 (^> c, Ao, zo) and the corresponding probabil¬ 
ity -P(0, c, Ao, zo) can be then calculated with the formulae 
similar to Eqs. (2) and (3), respectively. Therefore, the cal¬ 
ibrated luminosity is derived by 

-• _ f c f n L(fi,c,Ao,zo)P(n,c,Ao,zo)dfidc 
° ~ f c f n P(ti,c, A 0 ,z 0 )dtidc ’ (6) 


and its rms is estimated by 

2 _ fcfn&(Cl, c , A o,zo) - L( z o)] 2 P{n,c, A 0 ,z 0 )dCldc 
° ~ f c f n P(fi,c, A 0 ,z 0 )dtidc ’ U 


where L(zo) is the unweighted mean of Lo(f2, c, Ao, zo) in 
different c and Cl values. 


3 SIMULATION TESTS 

The current GRB samples that favor various luminosity indi¬ 
cators are very small, so that one cannot directly utilize our 
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approach to perform the calibration. The calibrations would 
nonetheless become possible in the future when enough data 
are accumulated. We therefore simulate a large sample of 
GRBs to examine our approach. The simulations aim to ad¬ 
dress the questions such as how many bursts are needed, and 
how narrow the redshift bin should be used, etc., given a par¬ 
ticular observed sample. We take the LZ-relation (Liang & 
Zhang 2005) as an example, which reads 

E^ ao = cE'“H'r, ( 8 ) 

where t b and E p are measured in the cosmic rest frame of 
the burst proper. We simulate 10 6 GRBs. Each simulated 
GRB is characterized by a set of parameters denoted by (z, 
E p , Ei ao , tb )• It is well known that the E v distribution of a 
bright BATSE GRBs presented by Preece et al. (2000) is well 
modelled by a Gaussian function. The HETE-2 and Swift 
observations of X-ray rich GRBs and X-ray flashes (XRFs; 
Heise et al. 2000; Lamb, Donaghy, & Graziani 2005) have 
considerably extend the E p distribution to a softer band. 
Liang & Dai (2004) studied the observed E p distribution of 
GRBs and XRFs, combined with both HETE-2 and BATSE 
observations, and found that the observed E p distribution 
for GRBs/XRFs is bimodal with peaks at ~ 30 keV and 
~ 200 keV. The ~ 30 keV peak has a sharp cutoff at the 
low energy end, likely being due to the instrument thresh¬ 
old limit. A recent study of a Swift burst sample marginally 
reveals such a bimodal distribution (Zhang et al. 2006). We 
therefore model the E v distribution by combining the obser¬ 
vations of BATSE and Swift, i.e. 


dp 

d log E p 


+ 


0.70 

0.56 yJV/2 
0.30 

0.56 yJV/2 


exp[— 2 ( lo g^- 2 - 30 
PL v 0.56 

^log E p — 1.55 
ex P[ 2 (- (TXr- 


( 9 ) 

( 10 ) 


with a cutoff at E p = 30 keV (see Fig. 1). The Ei so distri¬ 
bution is obtained from the current sample of GRBs with 
known redshifts. Since the Ei SO distribution suffers obser¬ 
vational bias at the low E lao end, we consider only those 
bursts with Ei ao > 10 51 ' 5 ergs, and get 2 dp/d\ogE lao oc 
—0.3 log E lao . The redshift distribution is assumed following 
the global star forming history of the universe. The model 
SF2 of Porciani & Madau (2001) is used. We truncate the 
redshift distribution at 10. A fluence threshold of S 7 = 10~‘ 
erg cnF 2 is adopted. 

We assume that these GRBs satisfy the LZ-relation and 
derive t b from the simulated E p and Ei ao . Since the observed 
tb is in the range of 0.4 ~ 6 days, we also require that tb is in 
the same range to account for the selection effect to measure 
an optical lightcurve break. Since the observed <r x /x is about 
10 % — 20%, the simulated errors of these observables are 
assigned as a x /x = 0.25 k with a lower limit of a x /x > 5%, 
where x is one of the observables E p , S 7 , and tb, and k is a 
random number between 0 ~ 1. Our simulation procedure 
is the same as that presented in Liang & Zhang (2005). 

With the simulated GRB sample we examine the plau¬ 
sibility of our calibration approach. We consider only a flat 


2 Our simulations do not sensitively depend on the /q so distribu¬ 
tion. We have used a random distribution between 10 51 ' 5 ~ 10 54,5 
ergs, and found that the characteristics of our simulated GRBs 
are not significantly changed. 


universe with a varying TLm- We picked up two samples with 
100 GRBs in each group 3 . The first group has a narrow red¬ 
shift bin (i.e. z = 2.0 ± 0.05) and the second group has a 
wide redshift bin (i.e. z = 2.0 ± 1.0). We then derive the 
parameters c, ai, and a 2 using the multivariable regression 
analysis (Liang & Zhang 2005) for different cosmological pa¬ 
rameter RIm) and evaluate the dependences of the derived 
parameters on SIm- The dependences of these quantities on 
Qm are quantified by the Spearman correlation, and the re¬ 
sults are presented in Figure 2. It is found that c strongly 
depends on SIm regardless of the value of Az, as is expected. 
On the other hand, while ai and 02 are strongly correlated 
with S2 for the case of Az = 1.0, they are essentially inde¬ 
pendent of SIm for Az = 0.05. These results suggest that 
once Az is small enough the influence of cosmological pa¬ 
rameters on both a 1 and 02 becomes significant lower than 
the observational uncertainty and the statistical fluctuation. 
This makes the calibration of both a 1 and (12 possible with 
a GRB sample within a narrow redshift bin. 

The selection of Az is essential to most optimally es¬ 
tablish the calibration sample. Two effects are needed to 
take into consideration to select Az, i.e. the observational 
errors of the sample and the statistical fluctuation effect. 
The most optimal calibration sample requires that the vari¬ 
ations of the standard-candle parameters caused by varying 
cosmology should be comparable to the variations caused by 
these two effects. In such a case we could establish a sample 
with a large enough N to reduce the fluctuation effect while 
in the mean time with a small enough Az so that the depen¬ 
dences of both a 1 and a? on Qm are not dominant. Since the 
relation between 01 (or CI 2 ) and log Qm is roughly fitted by 
a linear function (see Figure 2), we measure the dependence 
by the chance probability (P) of the Spearman correlation. 
If P < 10 -4 the dependence is statistical significant, and the 
sample is inappropriate for the calibration purpose. Figure 
3 shows the distributions of logP for a\ (left) and <12 (right) 
in the A z-N plane, assuming the current observational er¬ 
rors for the observables. The grey contours mark the areas 
that the dependences of ai and 02 on Qm are statistical 
significant. We find that P dramatically decreases as Az in¬ 
creases for a given N. Given a P value, Az initially decreases 
rapidly as N increases but flattens at N > 50. This indicates 
that the statistical fluctuation effect is much lower than the 
observational errors for a sample with N > 50. We can see 
that with the current observational precision, Az ~ 0.3 is 
robust enough to calibrate both ai and < 22 . Increasing the 
GRB sample size alone does not improve the calibration 
when N > 50, since the ai and 02 errors are dominated 
by the observational uncertainties in the data. In order to 
improve calibration further, higher observational precision 
of Ei ao , E p and tb is needed, which requires a broad-band 
7 -ray detector and good temporal coverage of the afterglow 
observations. 

The observed GRB redshift distribution ranges from 
0.0085 to 6.3. We examine if there exists a preferable red¬ 
shift range for the calibration purpose. We randomly select 
a sample of 25 GRBs at z c ± 0.3, and perform the multi- 
variable regression analysis to derived a\ and 02 from this 
sample by assuming a flat universe with Om = 0.28. The 


3 To avoid the statistical fluctuation effect we use a large sample. 
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Figure 1. The observed E p distribution: dashed line — derived 
from bright GRB sample (Preece et al. 2000); step-line — Swift 
data (Zhang et al. 2006); smoothed-curve — our model with bi- 
modal Gaussian distribution (Eq. 9). The dark region marks the 
cutoff at E p < 30 keV due to the instrument threshold limit. 

derived ai and a 2 are plotted as a function of 2 C in Figure 4. 
We find that they are not correlated with z c , and their vari¬ 
ations are essentially unchanged, i.e. ~ 0.15. This indicates 
that there is no evidence for a vantage redshift range to cal¬ 
ibrate ai and a 2 . It is therefore equivalent to select a sample 
at any redshift bin to calibrate a\ and ( 12 - Such a sample is 
likely to be established with GRBs at z c = (1 — 2.5), since 
the observed redshift distribution peaks in this range. The 
cosmological dependence is less significant at lower redshifts. 
So, a lower redshift (e.g. z c = 1) sample is preferred for c 
marginalization. 

According to Figure 3, the best strategy to perform 
GRB standard candle calibration is to establish a moder¬ 
ate GRB sample (e.g. 25 bursts) within a redshift bin of 
A z ~ 0.3 at a fiducial intermediate redshift (e.g. z c ~ 1 
or z c ~ 2). We simulate a sample of GRB with N = 25, 
Z = 1 ± 0.3, and derive ai and (22 as a function of f2 m 
in Figure 5. The calibrated a\ and 02 are 1.93 ± 0.07 and 
— 1.23A0.07, respectively, where the quoted errors are at 30 - 
significance level. 


4 CONCLUSIONS AND DISCUSSION 

We have explored in detail an approach to calibrate the 
GRB luminosity indicators, L(f2) = c(f2)Q(f2|X; A), based 
on the Bayesian theory without a low redshift GRB sam¬ 
ple. The essential points of our approach include, (1) cali¬ 
brate A with a sample of GRBs in a narrow redshift bin A z; 
and (2) marginalize the c value over a reasonable range of 
cosmological parameters for a given GRB sample. We take 
our newly discovered multi-variable GRB luminosity indica¬ 
tor Ei so = cEp 1 t 1 ^ 2 (LZ-relation) as an example to test the 
above approach through simulations. We show that while c 
strongly depends on cosmological parameters, both ai and 
0,2 do not if Az is small enough. The selection of Az depends 
on the size and the observational uncertainty of the sample. 
For the current observational precision, we find Az ~ 0.3 is 
adequate to perform the calibration. 


Figure 2. Comparison of the dependences of c, aj, and 02 on 
for a sample of 100 GRBs distributed in 2 = 2.0 ± 1.0 (left 
panels) and in 2 = 2.00 ± 0.05 (right panels), respectively. Cur¬ 
rent observational errors are introduced for the simulated bursts. 
The dependences are measured by the Spearman correlation, and 
the correlation coefficient (r) and its chance probability ( P) are 
marked in each panel. 
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Figure 3. Distribution of logP in the ( N , A 2 )-plane. The grey 
contours mark the areas where the dependences of ai and 02 on 
!1 m are statistical significant (P < 10“ 4 ). The white region is 
suitable for the calibration purpose. 


It is also found that the calibrations for both ai and 
02 are equivalent for samples at any redshift bin. The best 
strategy would be to collect GRBs within a narrow redshift 
bin around a fiducial intermediate redshift (e.g. z c ~ 1 or 
Zc ~ 2), since the observed GRB redshift distribution is 
found to peak in this range. Our simulation suggests that 
with the current observational precision of measuring GRB 
isotropic energy (Ei so ), spectral break energy ( E p ), and the 
optical temporal break time (tb), 25 GRBs within a redshift 
bin of Az ~ 0.30 would give fine calibrations to the LZ- 
relation. Inspecting the current GRB sample that satisfies 
the LZ-relation, we find that nine GRBs, i.e. 970828 (2 = 
0.9578), 980703 (2 = 0.966), 990705 (2 = 0.8424), 991216 
(2 = 1.02), 020405 (2 = 0.69), 020813 (2 = 1.25), 021211 
(2 = 1.006), 041006 (2 = 0.716), and 050408 (2 = 1.24) are 
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Figure 4. The variations of ai and a 2 as a function of z c . The 
calibration sample consists of 25 simulated GRBs at 2 C ± 0.3. 
£Im = 0.28 is adopted. 
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Figure 5. The variations of a\ and 0.2 as a function of The 
calibration sample consists of 25 simulated GRBs distributed in 
the redshift bin z = 1 d= 0.3. The dashed lines enclose the 3cr 
significance regions. 


roughly distributed in the redshift range z = 1.0 d= 0.3. We 
expect roughly 15 more bursts to form an adequate sample 
to calibrate the LZ-relation. 

The observed redshift distribution for the current long 
GRB sample covers from 0.0085 to 6.29. There have been 
suggestions that GRB properties may evolve with redshift 
(e.g. Lloyd-Ronning et al. 2002; Amati et al. 2002; Wei 
Sz Gao 2003; Graziani et al. 2004; Yonetoku et al. 2004). 
Among the proposed GRB luminosity indicators the cos¬ 
mological evolution effect has not been considered. With 
the current GRB sample with known redshifts, it is diffi¬ 
cult to access whether and how GRBs evolve with redshift. 
Nonetheless, since our calibration approach makes use of a 
GRB sample in a narrow redshift bin, the evolution effect es¬ 
sentially does not affect on the calibration of the parameter 
set A, the set of the power index (indices) in the luminos¬ 
ity indicators. However, it could significantly impact on the 
c-marginalization. 
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